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Abstract 

CN ! In this paper, we derive the quantum hydrodynamics models based on the moment 

closure of the Wigner equation. The moment expansion adopted is of the Grad type 
firstly proposed in [T5]. The Grad's moment method was originally developed for the 
' Boltzmann equation. In [7|, a regularization method for the Grad's moment system of 

the Boltzmann equation was proposed to achieve the globally hypcrbolicity so that the 
r~| ' local well-posedness of the moment system is attained. With the moment expansion 

of the Wigner function, the drift term in the Wigner equation has exactly the same 
moment representation as in the Boltzmann equation, thus the regularization in [7] 
applies. The moment expansion of the nonlocal Wigner potential term in the Wigner 
equation is turned to be a linear source term, which can only induce very mild growth 
of the solution. As the result, the local well-posedness of the regularized moment 
' system for the Wigner equation remains as for the Boltzmann equation. 

' Keywords: Moment closure; Wigner equation; Boltzmann equation; quantum hy- 

^SJ . drodynamics 

^ '■ 1 Introduction 



The Wigner equation was proposed by Wigner in 1932 as a counterpart in quantum 
mechanics of the Bohzmann equation [37j . Different from the distribution function in the 
Bohzmann equation, the Wigner function may take negative values, so it is not a true 
probability distribution function. But it is called a "quasi-probabihty distribution func- 
tion" in the phase space because it ahows one to express quantum mechanical averages 
(related to statistical moments) in a form which is very similar to that for classical aver- 
ages. Readers may refer to the review article by Hillery et al. [20] for more details of the 
properties of the Wigner function. It offers a convenient interpretation of the ensemble 
in the form of a quasi-probability distribution in the phase space defined by position x 
and momentum p and has been proved to of great use not only as a calculation tools 
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but can also provide insights into the connections between classical and quantum me- 
chanical mechanics. Actually, Lions et al. have given a rigorous mathematical proof that 
in the semiclassical limit, the Wigner equation gives the Vlasov (or Liouville) equation 
in the phase space [28]. As one of the most accurate equations in quantum mechanics, 
the Wigner equation has many advantages over the Schrodinger equation, the nonequilib- 
rium Green method and the density matrix method in simulating the carrier transport in 
semiconductor devices because the description of boundary conditions and collision oper- 
ator for the semi-classical Boltzmann equation can be extended to the Wigner equation 
due to their strong similarity [13] • There has been an increasing interest in the Wigner 
equation as the size of semiconductor devices goes into the nano-scale under which the 
quantum effects becomes not negligible any more [13]. In some devices, the quantum 
effect even plays a dominating role, e.g., the resonant tunneling diode (RTD) has been 
extensively studied recently using the numerical methods based on the Wigner equation 
[15j . The numerical method for the Wigner equations has attracted many researchers from 
different fields [3ll[I71[29l[25l[22l[26l|38l[5]. The deterministic numerical methods have 
been successfully used in simulating one-dimensional devices, but are not expected to be 
directly used for multi-dimensional devices simulation due to its formidable expense in 
memory storage and computational time. One practical approach to investigate a higher 
dimensional devices where quantum effects are relevant is to use quantum hydrodynamics 
models which are moment systems derived from the Wigner equation. Because the close 
connection of the Wigner equation and the Boltzmann equation, many moment meth- 
ods devised for the Boltzmann equation have been extended to the Wigner equation, see 
[16t [T2] and references therein. Equations derived from the Wigner equation are called 
quantum drift-diffusion equations, quantum Euler equations and quantum hydrodynam- 
ics equation, and numerical simulations based on such moment equations are extensively 
studied [39l IHl 1211 [H]- In this paper, we will extended the moment method recently 
proposed in [7] for the Boltzmann equation to the Wigner equation. 

In 1940s, Grad [18] proposed a moment method to approximate the Boltzmann equa- 
tion, and a 13-moment model is given as an extension of the classic Euler equations. 
However, this model was quickly found to be problematic [1^. Its major deficiencies in- 
clude the appearance of subshocks in the structure of a strong shock wave and the loss 
of global hyperbolicity. Later on, a number of regularizations were raised to solve or 
alleviate these problems. Levermore [27] introduced a promising way to achieve global 
hyperbolicity, while 14 moments are needed and the explicit expressions of the equations 
cannot be written. Jin and Slemrod [23] gave a regularization of the Burnett equations via 
relaxation, which resulted in a set of equations containing the same variables with Grad's 
13-moment theory, and no subshocks appear in the structure of shock waves. Struchtrup 
and Torrilhon [33] regularized Grad's system by integrating the moment method with 
Chapman-Enskog expansion, and the resulting system is called as the R13 equations. 
And in |32j . Struchtrup improved the R13 equations by the "order of magnitude" method. 
The R13 system also removes the discontinuities in the shock wave, and it extends the 
region of hyperbolicity considerably [M] ■ Recently, Torrilhon [5S] tackled the problem of 
hyperbolicity by introducing the Pearson- Type-IV distributions, and the resulting system 
is proven to be able to gain a much larger hyperbolicity region. 

Due to the complexity of the explicit expressions, systems with large number of mo- 
ments are not investigated until recently. In [36], Torrilhon and his coworkers developed 
a software named ETxx [2] which is able to generate moment systems with almost any 
number of moments in one-dimensional space. And in [3], some numerical results for a 
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shock tube were carried out to show the behavior of characteristic waves in the extended 
thermodynamics. A numerical method solving large moment systems was proposed in [8], 
and therein the regularization technique in [33] was applied to general moment systems. 
In [To], the order of magnitude method was also integrated into large moment systems. 
In [6], the authors focused on the one-dimensional velocity space, and the characteris- 
tic polynomial of the quasi-linear coefficient matrix is found to be very simple; thus a 
brand-new regularized model with global hyperbolicity is proposed by the correction of 
the characteristic speed. Such regularization is extended to the multi-dimensional velocity 
space in [7]. 

It is clear that the moment expansion of the drift term for Boltzmann equation can be 
extended to the Wigner equation. The resulting convective terms in the moment system 
expanded from the drift term of the Wigner equation has exactly the same format as that of 
the Boltzmann equation. Thus the method of the hyperbolic regularization in [7] is applied 
to the Wigner equation, too, to achieve the global hyperbolicity. As the difference of the 
Wigner equation from the Boltzmann equation, the nonlocal Wigner potential term due to 
the electric potential is also expanded using Hermite polynomial. It is interesting for us to 
find that this term can be represented in the moment expansion style with very compact 
expressions. The overall formation of the nonlocal Wigner potential term in the moments 
is turned into a linear source term, with compact sparse coefficient matrix. Moreover, 
the contribution from the lower order moments in this term is always on the higher order 
moments. Immediately, the coefficient matrix in the linear source term is strictly lower 
triangular, thus it is a nilponent matrix. Noticing that the moment expansion of the 
relaxation scattering term produces a linear source term providing an exponential decay 
of the high order moments. As a result, the growth of the high order moments in time due 
to the source term is essentially slower than exponential growth rate. This makes that the 
derived moment system is formulated as a quasi-linear system, plus a linear source term 
which can induce only very mild growth of the high order moments. Since the convection 
term in the system is guaranteed to be globally hyperbolic by the regularization, the local 
well-posedness of the system is out of box. 

The rest of this paper is arranged as follows: in Section 2 we present the elementary 
formula of Wigner equation. The moment expansion of Wigner equation is carried out in 
Section 3 and the system obtained is closed by truncation of the expansion and regularized 
using method in [7] in Section 4 to achieve the final hyperbolic system. In Section 5, we 
discuss the simple case in ID space for better understanding of the structure of the derived 
moment system. Concluding remarks are in the last section. 



2 Wigner Equation with Smooth Potential 

We consider a one-particle system in a statistical mixture with N states described by 
the wave functions = x), each with a probability ^ 0, with k = 1, - ■ ■ , A^statej 
which satisfies Efc=T ' Pk = I- The wave functions "^k^t, x) of state k obeys the Schrodinger 
equation 

where i = h is the reduced Planck constant, m is the particle effective mass assumed 

to a constant in this paper, and V{t,x) is the electric potential energy (which will be 
called potential for short hereafter). We can also regard the one-particle system as a 
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many-particle system with interpreted as the percentage of the particles occupying 
state k. This is the general way the Wigner equation used in modeling carrier transport 
in semiconductor devices. We construct the Wigner function f{t,x,p) which is a quasi- 
probability density function in the phase space {x,p), as usual, by 

-| -Estate /» 

f{t,x,p) = —— V Pk / nit,x + y/2)^kit,x-y/2)e'yp/'^^d^y. (2) 

Using the Schrodinger equation ([T]) and the definition of the Wigner function given above, 
we can derive the Wigner equation (or the quantum Vlasov equation) as been done in 
[S71I2D] 

^ + --V,/ + e[y]/ = 0, x€R',p€R^ (3) 
at m 

where the nonlocal Wigner potential term 0[V]/ is a pseudo-differential operator. It is 
defined by 

{e[V]m,x,p)= [ VUt,x,p')fit,x,p-p')dp', (4) 



where the Wigner potential is as 

Vu){t,X,p) 



{2-Khfh 



For convenience, we simply take m = 1 later on. Similar to a Vlasov-Poisson system, a 
Wigner-Poisson system could be also considered to include the self-consistent electric field 
induced by the redistribution of the electrons. In this paper, we focus on the case that 
V{t^x) is a given smooth function. The pseudo-operator can be written into an 

equivalent form [20] 

(em/)(t,.,p) = -x;^^^^^, (5) 

where A = (Ai,A2,A3), |A| = ELi ^i, = nii^;^, 



rr ^ _ rr 



dx^ fi^rr^' dp^ 11 

and the summation over A has to be extended over all non-negative integer values of 
Ai,A2,A3 for which |A| is odd. In the semi-classical limit /i — )■ 0, 0[y]/ converges to 
the usual operator —VxV • ^pf- The equation ([Sj) describes particles movement without 
scattering corresponding to ballistic transport. The scattering effect may be considered 
by adding a scattering term to the right hand side of (l3|). So one obtains the following 
Wigner equation with a scattering term: 



, £c e R^, p G M^. (6) 

Scat 



The time-relaxation approximation of the scattering term is often used for the Wigner 
equation [15], and it has the same form of the BGK scattering term used for the collision 
of gas [5. It is easy to be study analytically and takes of the following simple form 

^/ _ /eg ~ / /_s 
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where r is the relaxation time and /eq is the equihbrium distribution. For example, 
when the density of electrons are not extremely high, we can assume that the equilibrium 
distribution is a Maxwellian distribution, 

where p{t, x) is the number density of particles at position a;, ks is the Boltzmann constant, 
T(t, x) is the particle temperature, and u{t, x) is the average momentum of particles. 
These macroscopic variables are related with the distribution function as below: 



p{t,x)= / f{t,x,p)dp, (9) 

p{t,x)u{t,x) = pf{t,x,p)dp, (10) 

Jm.3 

p{t,x)kBT{t,x) = [ \p-u\'^fit,x,p)dp. (11) 

The conservation of mass, momentum and total energy are all valid for the Wigner 
equation. Multiplying the equation ([6]) by 1 and p, direct integration with p and x gives 
us 

4/" fit,x,p)dxdp = 0, tGR+, (12) 

4- f pf{t,x,p)dxdp = - [ p{t,x)^^^^^dx, teR+. (13) 
dt Je3x]r3 J ^3 ax 

Multiplying the equation ^ by |pp and integrating by parts, we get the conservation of 
the total energy for the system Q and ([5]): 



/ f{t,x,p)\p\'^dxdp = -2[ p{t,x)u{t,x) ■VxV{t,x)dx, t£ 
drjM3xM3 Jr3 



(14) 



3 Grad Moment System 

In this section, we derive the moment system of the Wigner equation using the Grad 
type moment expansion. 

3.1 Hermite expansion of the distribution function 

Following the method in [SJ [TU], we expand the distribution function into Hermite 
series as 



f{t,X,p) = J2 Ut,^)'HT,a f^7==) 



(15) 



where a = {ai,a2,cx3) is a three-dimensional multi-index. The basis functions 'HT,a are 
the 3-dimensional Hermite functions defined by 

nrA^) = Jl-^T-^He^,iCd)ew ("f ) , (16) 
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where Hen{x) is the Hermite polynomial of order n 

He^ix) = (-l)" exp (^) ^ exp 1^-^^ . (17) 

For convenience, Hen{x) is taken as zero if n < 0, thus T-LT,a{$,) is zero when any component 
of a is negative. The parameter T in the expansion is the scaled local temperature as 

T{t,x) = kBT{t,x). (18) 

It is clear that the equilibrium distribution /cq is coincidently equal to the first term 
of expansion, i.e., 

Uit,x,p) = fo{t,x)nr,o ( P-gl^) ] , (19) 

where /o(i, a^) = /o(i, a^). So in this case the relaxation time approximation scattering term 
(j20p can be written as the linear combination of the basis functions of orders greater than 
or equal to 1, 

The definition of the Hermite function ()16p shows that each basis function is an expo- 
nentially decaying function multiplied by a multi-dimensional Hermite polynomial shifted 
by the local macroscopic momentum u and scaled by the square root of the local temper- 
ature T. 

If one uses arbitrary known function ii'(t, a;) and ^'(t, x) in (|15p to expand the distri- 
bution function f{t,x,p) as 



f{t,x,p) = at,x)nr',a (^7==-) 



(21) 



then the following relations between the macroscopic quantities u, 7~, the coefficients fa, 
and u', T', the coefficients can be derived as follows. 



P = fo = /o, (22a) 
p« = p«' + (4)J=i_2,3, (22b) 

3 

p\u - u'\^ + 3pT = Y.{T'f'^ + 2/^,J, (22c) 

d=l 

where is the unit vector with its d-th. entry to be 1. It is clear that the coefficients fa 
expanded using parameters u and T satisfy the following conditions: 

3 

/e. =0, J^/2e, = 0, i = l,2,3. (23) 

d=l 

Moreover, if we define the heat fiux qi and the pressure tensor P = {pij}, i,j = 1,2,3 
with 

Qi = l I \p-u\^{Pi-Ui)fdp, (24) 
2 Jm.3 

Pij = {Pi- Ui)iPj - Uj)fdp, (25) 
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then direct calculations give us the relations between them and the coefficients fa in (|15p 
as 
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(26) 



d=l 



Pij - ^^ij ^Pdd = (1 + 6ij)fe,+ey 



(27) 



By the definition of the temperature (jlip and (jlSp and the definition of the tensor pressure 
(j27|) , the scaled temperature T is a linear combination of pij as 



3.2 Moment expansion of the Wigner equation 

Now we are ready to derive the moment system by taking the moments of the Wigner 
equation. The general method to get the moment system is to first multiply the Wigner 
equation ^ of by polynomials of momentum p of different order and then integrate both 
sides over momentum p on M'^. One equivalent way is as follows. First, we substitute 
the expansion of the Wigner function ()15p into the Wigner equation ([6]), then we collect 
the coefficients of basis functions of the same order at both sides, and finally we equate 
the coefficients of the basis functions of the same order on both sides to yield the derived 
moment system. We plug the Hermite series (|15p into the Wigner equation ([H]), and make 
calculations by noting that the Hermite function ()16p used in this paper depends also on 
the time t and position x through u{t,x) and T{t,x), which is different from the general 
expansion using the Hermite functions depending only on the momentum p [30]. For 
convenience, we list some useful relations of Hermite polynomials as below [1]: 




(28) 



d=l 



With the relation (j27p . we then have 



Pij = SijpT + (1 + 6ij)fe,+ey 



(29) 




2. Recursion relation: Hcn+iix) = xHen{x) — nHen-i{x); 

3. Differential relation: He'^{x) = nHen-i{x). 



And the following equality can be derived from the last two relations: 



[Fe„(x)exp(-xV2)]' 



Hen+i{x) exp(-x^/2). 



(30) 



Especially, we have 




(31) 



With these relations, the part 




7 



of is expanded as 



r 



d=l 



(1=1 



+ 1^ (T/a-ed-e,- + ^^j/a-e, + (^j + l)/a-e,+e,) 



1 9r 



p — U 



T 



(32) 



Then using pip, we calculate the pseudo-operator term expressed in ([5]), and obtain 



ie[v]m,x,p) = J2 



a, A 



A! 'dx^ 



p — U 



(33) 



where the summation over A means the same as that in ([5]), and fa-x is taken as zero 

df 

when any component of a — A is negative. Finally, the scattering term — given in 

at 



becomes 



Scat 



la|>2 V V / / 



(34) 



noticing that /e^ = 0, d = 1, 2, 3. 

Collecting the three terms ([32]) . ([33]) and (f3l]l yielded after the substitution of the 
Hermite expansion p5|) into the Wigner equation ([B]), we can get the following general 
moment equations with a slight rearrangement by matching the coefficients of the same 
weight function: 



dfa 



•JJa . sr^ I dud uud 1 L I u I s-^ '^'lY^^ 



1 / dT 



dT 



dt 



d=i 



dt 

3 



(9x 



2 5t 



d=l 



+ ("i + l)/a-ed+ei) 



1 



5/o 



^mi I 9W ^ (V2i)l^l"i OV^ 
-H(|a| - 2)/^ - 2^ T^/a-A, 

A 



A! dx^ 

where H(x) is the Heaviside function defined by 

H(x) = 



0, if X < 0, 

1, if X ^ 0. 



(35) 



(36) 
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By setting a = in ()35p . we deduce the mass conservation 



By setting a = e^, with d = 1, 2, 3 and noting that /e^ = in (|35|) . we obtain 

I dud ^ dud\ dT dp ^ dfe^+e, dV 

which is simphfied as 



dud A dud\ dpjd dV 



P 



By setting a = 2ed, with (i = 1, 2, 3 and noting that /e^ = 0, we obtain 

+ E -.-^ + (1 + = 



(40) 



Noting that Yld=i ~ 0) we sum the upper equations over d to get 



3 



dT dT \ 2 v-^ / dqj Sn^ \ „ / . . n 

7^7 + E-.7^r + E 7t + T.P^^7^ = 0- (41) 



Since pT = ^ J2d=iPdd, we have 

— = — V - -— j = 1 2 3 (42) 



Substituting p9|) . (|4T]) and (j42|) into (|35|) . we ehminate the time derivatives of u and T 
and the spatial derivatives of T. Then the quasi-hnear form of the moment system reads: 

3 



^ + E(^^+-.9^ + (". + i)^J 



3 3^^/ 3 

+ E E T/a-e^-e, + (ttj + l)/„-e^+e,. " E ■^"-2^ 

j=l d=l ^ \ ^ k=l 



3 3/* ^ / 3 \ 3 4^ 



+ E (^-^4 + ^ 15 ^ J 15 ^^^--'^^-^^ + ("^- + 1)^-2^^+^.) j 

= --H(|a| - 2)/„ - g^fa-x, V|a|>2. 
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We should point out again the summation over A is extended over all the non-negative 
integers A^, d= 1, 2, 3 for which | A| is odd and greater than 1. Especially for |a| < 2, the 
moment system (|43p derived from the Wigner equation is the same as that derived from 
the Boltzmann equation [9]. 

With ([29|) , we can have the equations for pij by ([I3]l . Precisely, we have the equation 
for pii/2, i = 1,2,3, as 

i=i ^ i=i ^ ^ J j=i d=i ■' 

j=l \ d=l / 

If i 7^ j, we have pij = fa+ej, thus its equation is already in (jl3|) . 

We collect the equations (f37|) , ([39]) , (jl^ and (03]) together to obtain a moment system 
with infinite number of equations. Noting that the relation between Ud and fe^ given in 
()22p and the definition of qi and of pij given in ()27p . we can see that what we obtain is 
a quasi- linear system for fa- We would like to point out that the only difference between 
the system derived from the Wigner equation and that from the Boltzmann equation is 
the term with high-order derivatives of the potential V{t, x), which is a source term of the 
quasi- linear system of fa- 



4 Moment Closure with Global Hyperbolicity 

The moment system derived from the Wigner equation consists of (|37p , ()39p , (|44p and 
()43p . It is clear that this is a system with infinite number of equations taken p, Ud, Pij 
and fa, |a| ^ 3, as unknowns. To obtain a system with finite unknowns, we will truncate 
the expansion (jlSp and close the system following the method in [Tj. 

With a truncation of (|15p . (|43p will result in a finite moment system. Precisely, we 
let M ^ 3 be a positive integer and only the coefficients in the set Ai = {fa}\a\^M 

considered. Let Fm{u, T) denotes the linear space spanned by all 'HT,a {^^^^^^ with 
I a I ^ M, and the expansion (|15p is truncated as 

f{t,x,p) « Ut,x)nT.a , (45) 

with f{t,x,p) G Fm{u,T) and fa & -M- The moment equations which contain dfa/dt 
with \a\ > M are disregarded in (jl3]). Then, ([37]), (jM]) and (03]) with 2 < |a| < M 

lead to a system with finite number of equations. 
Following [7], we let 

Sm = {a G I \a\ < M}- 

Then for any a € 5a/, let 

i=i ^ * 
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to be the ordinal number of a in 5m, and the cardinal number of set Sm is 

which is total number of moments if a truncation with \a\ < M is applied. 
Let w = {wi, • • • , € and for each i, j G {1, 2, 3} and i / j, 

wi = p, W(eO = Ui, (47a) 

W(2e,) = W(e,+e,) = Pij, (47b) 

= fa, 3 < |a| < M. (47c) 
The moment system ([37|) . ([39]) . (fli|) and (|l3]l is collected in quasi- linear format as 

j=l J 

by taking the derivatives of fa+ej, \ct\ = M to be zero, where Mj and G are N x N 
matrices. The entries of Mj are given as the coefficients of the terms in (|37p . ()39p . ()44p 
and (j43p with derivatives of w. The entries of G arise from the nonlocal Wigner potential 
term and the scattering term. From (|44p. one can observe that 

GAr(2e,),Af{2e,) = -^{^ij -^), i,J = 1,2,3. (49) 

And (|43p indicates that the diagonal entries of the lower right part of G are 

Gj\f(a),j\fia) = H(|a| — 2), for \a\ > 2 and a / 2ej, i = 1,2,3. (50) 

From (IHUI). it is clear that 



1 dV 

G^(e.),i = --5^' ^ = 1'2,3. (51) 
The other nonzero entries of G from the nonlocal Wigner potential are as 

G^(„),Ar(a-A) -Q^, (52) 

where |A| is odd and |A| is greater than 1 and |a — A| 7^ 1 or 2. In case of |a — A| = 1, 
we have G_/\^(q) _/\/(q„;^) = since fe^ = 0, i = 1,2,3. In case of |a — A| = 2, we have 
(f52]) if a — \ = ei + Cj , i ^ j , since wj^(^ei+ej) = Pij = fei+Sj- The difference is in case of 
a — A = 2ei, i = 1,2, 3. Precisely by ([221), 'we have 

GMia)M2e.) --2^[^v-^) («_2e,)! 9^-2^^ ^^3^ 



All other entries of G vanishes except for the ones specified above. It is clear that 
the non-zero entries defined by (|49p and ()50p by the scattering term provide us a linear 
damping of the high order moment. The contribution of the nonlocal Wigner potential 
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matrix G is strictly lower triangular. Thus these part of the matrix is nilponent. As a 
result, the growth of the high order moments in time due to the nonlocal Wigner potential 
term is essentially slower than exponential growth rate. 

In ipHj) . we following Grad [18] take dfa+ejdxj, \a\ = M, as zero to make the system 
to be closed. It has been pointed out in [7] that it is not appropriate to set dfa+cj /dxj = 0, 
\a\ = M, as the closure proposed in |18j since the system is lack of hyperbolicity if the 
distribution function is far away from the equilibrium. To obtain a system with global 
hyperbolicity, we have to adopt the regularization given in [7]. For any a with \a\ = M, 
we define 



and 



n{j{a) = (aj + 1) 



dw , , dw 



^— 1 -I 




(54) 



Mj— — ='M.j— Tl\j{a)Ij^(^a)i for any admissible w, (55) 

■' \a\=M 

where Ik is the A;-th column of the N x N identity matrix. We regularize the system ()48p 
as 

j = l J 

which is the quantum hydrodynamics model we derived. It has been proved in [7] that 



Theorem 1. The regularized moment system (|56p is hyperbolic for any w with positive 
temperature. Precisely, for a given unit vector n = (ni,n2,n3), the matrix 



3 

E 



nj'Mj{w) (57) 



is diagonalizable with eigenvalues as 

u-n + Cn,mVT, l^n^m^Af + 1, (58) 



where Cn,m is a root of m-order Hermite polynomial, and satisfies Ci^m < ■ ■ ■ < C 
The structure of the N eigenvectors can be fully clarified. 



m,m ■ 



Based on this theorem, the regularized moment system (j56p is locally well-posed due 
to the hyperbolicity. We would like to mention here that the regularization here actually 
does not add any new terms to the system ()48p . On the contrary, it has erased the terms 
in ()43p with a factor aj + 1 in its coefficient for the equations of fa with |a| = M only. 



5 Regularized Moment System in ID 

In ID space, the structure of the system we derived is significantly simpler than in 
multi-dimensional case. In this section, we give the detailed formation of the moment 
system in ID case since the ID Wigner equation is already very useful in modelling of 
RTDs US]. 



12 



The ID Wigner equation is as 



A=l,3,5,--- 

where effective mass is again assumed to be 1 for convinence. The equilibrium distribution 
feq is assumed to be a ID MaxwelHan distribution 

The ID Wigner function is expanded as 

f{t,X,p) = ^fn{t,X,p)'HT,ni^^-7=^], (61) 
neN ^ ^' ' 

where 

UT,n{x) = ^r-("+i)/2/7e„(a;) exp(-xV2), (62) 
V 27r 



and the ID density p and average momentum u are defined by 

/oo roo 
f{t,x,p)dp, u = - pf{t,x,p)dp, (63) 
-oo P J —oo 

and the scaled temperature T defined by 

1 1'°^ 

T=kBT=- {p-uffit,x,p)dp. (64) 

P J —oo 

The above equation plus the definition for the pressure ()25p for P shows that 

P = pT. (65) 

In one dimension, (|23|) becomes 



/i = 0, /2 = 0. (66) 

The relation between the heat flux q and the coefficient /s of the expansion is turned 
to be 

9 = 3/3. (67) 

We present below the details of the regularized moment system truncated at very low 
order based on (|56|) . though the closed moment system up to any order can be written 
in a general formation. We will demonstrate that with the low order expansions of the 
distribution function, the moment system obtained is able to capture very typical quantum 
effects. As the simplest closed regularized moment system truncated up to M = 3 

is 

dp dp du , , 

du du dP dV 

'm^p''d-x^^ = -^P^ 

dP/2 ^dP/2 3pdu 3^/3 _Q (j^. 
dt dx 2 dx dx 
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+ = _i/3 + 

dt 2p'^ dx 2p dx dx " ^ r,. o •^Pi 



r 



24 9x3 ^ 



(71) 



^2 d^v 

where tt^P is the quantum correction term yielded by the nonlocal Wigner potential. 

24 ox'^ 

We reformulate (f68|) . ([69]) . (fTOj) . ([TT]) into a matrix form 



where w = {p, u, P/2, f^)^ and 



M 



dw , ^ , , dw 



/ u 





Gw. 



(72) 



V"2p2 



G 





avi 

p 










o\ 




2 





u 






P 




3P 


u 


3 





P 










P 










\ 
























1 








-rl 



(73) 



(74) 



^24 9x3 

Here we consider the following example to show the difference between the quantum 
moment system and the classical moment system. The Boltzmann equation in ID with a 
potential 

exp 



V{x) 



0, 



if |x| 
else, 



< 1, 



admits a steady solution 



f{t,x,p) 



exp 



with its moments as 



p = P = exp(-y(x)), 



/3 = 0. 



(75) 



(76) 



(77) 



p, u, P and /s in ()77p satisfy the moment system derived from the Boltzmann equation 

which can be obtained from (I72p by removing the term tt^P- The solutions given in 

24 ox-^ 

(j77p show that the classical behavior of particles obeys the Newtons law, i.e., the particles 
with low kinetic energy (p^/2 at position x ) can not move to a new position with energy 
V{x') which is greater than + V{x). As a result, p in (j77p will not change in time. 

In the quantum mechanics, the particles behave like a wave that the particles with 
lower energy can be tunneled through a potential barrier. The moment system derived 
from the Wigner equation will reflect this point. Let us consider the system ()72p with 
([77P as initial value. We carry out an asymptotic analysis to study the behavior of the 
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solution around t = 0, using time t as the asymptotic expansion parameter. The ansatz 
we adopted is as 



/3 = /3,0 + /3,it + O(i'), (78) 

P = Po + Pit + P2t^ + 0{t^), (79) 

U = Uo + Uit + U2t'^ +U3t^ + 0{t^), (80) 

P = PO + Plt + P2t^ + P3t^ + P4t^ + O(t^), (81) 

where the leading order terms are as the initial value 

po = Po = exp(-y(x)), no = /g.o = 0. (82) 
Substituting the ansatz into ()7ip . we immediately have 

/3,i = ^^('Hx)po, (83) 

thus /s is as 

h = Y^V^\x)p,t + 0{f) (84) 
Plugging dHlD and (179D - dill into ([70]), we obtain 

A = 0, ^2 = -|^(f(=^)(x)po), (85) 



which yields 

P = Pn - 

8 dx 



P = Po - f A ( F(3)(:,)po) + 0(i=^). (86) 



Similar procedure gives us 

d2 



U = Uq + 



24 dx2 V 



(y^^\x)p,)t' + 0{t^), (87) 



d3 

P = Po 



Noting that po = exp(— y(x)), we can rewrite ([88|1 into 



p = P^ + ^g{x) + 0{t^), (89) 

where 

ff(x) = -^(l^(^)exp(-2F)), (90) 
whose figure is plotted in Figure [TJ 

It? d '^V 

From ([8^ . we see that the heat flux /a is changed by the term — po, which is the 
difference between the Boltzmann equation and the Wigner equation. The change of heat 



flux induces the change of the temperature (j86p thus the electrons in part of the domain 
become cooler than the electrons in the else part of the domain. This makes the electrons 
with higher temperature moving to the domain with lower temperature ()87p . eventually 
resulting the redistribution of the density of the electrons (|88p . (|89p shows that how the 
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show only g(x) > 

-0.5 0.5 

J I I I I I I I I I I I I I L_ 



le+05 - 



le-05 



le-10 




- le+05 



le-05 



le-10 



Figure 1: The figure of g{x) in ()90p . The curve of g{x) is plotted in semi-log scale only in 
the subintervals where g{x) > 0. g{x) = for x is outside of (—1, 1) has not been plotted 
for it is obvious that g{x) = 0. 



density p changes at a very small t, and whether it will increase or decrease depends on 
the sign of g{x), which is shown in Figure [TJ So it is observed from Figure [1] that the 
particles will be redistributed because of the quantum term. Precisely, at a very small 
time t, the electron density in the subintervals of (—1,1) where g{x) > 0, whose curves 
are plotted in Figure [U is increasing, and the electron density in the rest subintervals of 
(— 1, 1) where g{x) < without its curves plotted, is decreasing. 

In the remain part of this section, we present the ID system obtained by truncating 
the expansion ()6ip at M = 4, 5 and 6 for the reader's convenience to carry out numerical 
simulations. If we truncate up to M = 4, we obtain a closed system consisting of ([68|) - 
roD and 

-P, 



dt dx 2 dx dx dx 



h , 

' r 24 9x3 



and 



dt p dx dx dx T 



(91) 



(92) 



And the system is collected into the matrix form (|72p with w = (p, n, P/2, f-^, f^)'^ and 



M 



/ 


u 


P 








o\ 







u 


2 














P 











3P 


u 


3 







p2 


4/3 


P 












u 


4 




2p2 













1 


P 




\ 





P 


P 


7 



(93) 
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\ 


dVl 

dx p 





























24 dx^ 








1 
r 




1 















\ 

Similarly we can write out the system truncated at M 
(UnD, and 



(94) 



5 which consists of 



and 



dt 



at 



hdP du 
p ax ox 



2 dx dx dx dx 



hdP_ 
p dx 



p dx dx dx 



d^V 



h 

' T 1920 dx^ 



T ' 



fo, 



which is expressed into the matrix form {[72 

/ 



M 



u 



2«2 

' 2p2 




P 
u 

3P 

4/3 
5/4 





with where w = (p, u, P/2, /s, h, /s) 




2 

P 

u 

P 

4 
"a 



T 



(95) 

(96) 
and 









o\ 











3 








u 


4 





P 








u 


5 








k 


P 


u] 


p 


P 



(97) 





dV l 

dx p 


d'^V 
24 5x3 


h'^ d^V 
V~ 1920 dx^ 

In the case of truncating at M = 6, the moment system obtained consists of 
dZni), dSI]), (USD and 















\ 




































1 
























r 


1 
































r 












1 














-rl 



(98) 



dt 



dh hdP ^hdh ^^.du dT 

^ S 1" "^5"n~ + 3/4^ 

dt p dx p dx dx ox 

, T-9/4 , dh , „9/6 h 

ox ox ox T 

hdP_3hdh^ 1^7-/3 + + u^ = 
p dx p dx 2 dx dx dx 



d^V 
1920 3^5" 



/oi 



h . h'd'v 

r 24 dx^ 



(99) 



(100) 
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This system can be formulated in the matrix form of ()72p with w 
and 



M 



u 



3P/4 



\ 2p^ 



P 

u 

3P 
4/3 
5/4 
6/5 




V 





dVl 

dx p 


24 dx'^ 


a^y 

1920 9x5" 




Ph 
2p2 












2 

P 
u 

P 

4 

2p 
P 



p 




1 

r 








3 

u 
P 

p 

p 





1 

T 
















4 

u 5 
P 

— u 
P 





{p,u,p/2j3,h,h,hy 

o\ 












1 

T 





p 

p 



7 
\ 








(101) 



(102) 



24 9x3 

Observing M's in ([12]), ([SZD, (fTUTT) and G's in (HH), ([MD, &, (HMD, we find that 
the left hand side of the moment system of ()72p is the exactly the same as the moment 
system derived from the classical Boltzmann equation, and the quantum mechanical terms 
with a typical characteristic which involves the high order derivatives of the potential V 
are only appear in the matrices G in the right-hand side of (j72p . It is clear that every 
G is lower triangular. The formation of matrices G shows that the quantum potential 
term works in the way by letting high-order moments get information from the moments 
of different lower orders, as low as to order 0. 



6 Conclusion 

We extend the moment closure method |7j for the Boltzmann equation to its quantum 
counterpart, the Wigner equation. And we obtain a hyperbolic moment system, i.e., 
quantum hydrodynamic system. The method to derive the quantum hydrodynamic system 
with the global hyperbolicity is systematic that systems for arbitrary number of moments 
are obtained at once. Numerical simulations are to be carried out to demonstrate that 
the quantum effects are able to be captured by using the quantum hydrodynamic model 
derived. 
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